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In presence of a small magnetic field h, the elementary excitations in the scaling two-dimensional 
Ising model are studied perturbatively in h in the ferromagnetic phase. For excitations with large 
numbers n, the mass spectrum is obtained in the first order in h. The decay widths of excitations 
with energies above the stability threshold are calculated in the leading /i 3 -order. 
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During the past decades there has been much progress 
in study of the Ising Field Theory (IFT). For a recent 
review, see Q. Being the field-theoretical realization of 
the two-dimensional Ising model in the scaling limit, IFT 
provides a direct information about the 2D Ising univer- 
sality class. IFT can describe also dynamical properties 
of many one-dimensional condense-matter systems, the 
best known example is the one-dimensional Ising magnet 
at T = near the quantum phase transition point 0, Q • 
Viewed as a particle theory, IFT turns into a beauti- 
ful model of quark confinement in which new ideas and 
methods can be efficiently worked out 0, 0, 0, ■ 

Physics of IFT is determined by the scaling combina- 
tion 77 — m/|/i| 8 / 15 of parameters m and h, which are 
proportional to the deviations of temperature and mag- 
netic field in the 2D Ising model from their critical point 
values T = T c and H = 0. Being not integrable for 
generic 77, the IFT has exact solutions only at 77 = ±00 
and 77 = 0. The former case is free-fermionic and corre- 
sponds to the Onsager's exact solution Q of the lattice 
Ising model in zero field. Solution of IFT at 77 = 0, that 
describes the scaling Ising model at the critical tempera- 
ture in a magnetic field, was found by Zamolodchikov || . 
Beyond the integrable cases, the main progress has been 
achieved by exploiting perturbation expansions around 
77 = and 77 = ±00 H H H [H [H E3- Irithe present 
paper we focus on the perturbative study of the Ising 
field theory in the limit 77 — ► +00, i.e. at m > and 
h — ► 0. This case corresponds to the ferromagnetic (low 
temperature) small-field scaling regime of the 2D Ising 
model. 

At h = 0, m > 0, the particle sector of the IFT con- 
tains one free fermion of the mass m, which can be inter- 
preted as a domain wall (kink) interpolating between the 
two degenerate vacua. External magnetic field h > re- 
moves degeneration, and induces a long-range confining 
interaction between fermions (kinks). As the result, an 
isolated fermion gains an infinite energy at a whenever 
small field h, and the kink-antikink bound states become 
the single-particle excitations of the model. This makes 
a link with the quark confinement physics. Accordingly, 
in what follows the fermions and their bound states are 
referred to as " quarks" and " mesons" , respectively. At a 



small external field, the weak confinement regime is real- 
ized: the meson masses M„ fill dense the interval above 
2m, the mesons with M„ > 2Mi are unstable, where Mi 
is the lightest meson mass, Mi « 2m. 

Though the qualitativepicture outlined above was un- 
derstood already in 1978 only few quantitative results 
on the weak confinement in IFT have been obtained ever 
since. The initial part of the meson spectrum M„ in the 
region 



M n — 2m <C m 



(1) 



was first obtained by McCoy and Wu jj| in the non- 
relativistic approximation. Relativistic corrections to 
this spectrum in the same region Q have been recently 
found by Fonseca and Zamolodchikov [(J ^J- In the 
present paper we calculate perturbativly in h the meson 
masses M„ for large n » 1, as well as the decay widths 
of the nonstable mesons. 

Ising field theory is determined by the euclidean action 
describing massive interacting Majorana fermions in the 
plane K 2 . For its explicit expression see for example Refs. 
[lj, |a] . The corresponding hamiltonian can be written in 
the form 



n 

where V 



dp 
2^ 



w(p) 4 a p + v > 



dxcr(x), 



(2) 



(3) 



and u(p) = \/p 2 +m 2 is the spectrum of free fermions. 
Fermionic operators a£, a p creating and annihilating a 
bare quark with a momentum p obey the canonical anti- 
commutational relations 



{ a P, %>} = 27 r5(p-p'), {a p ,a p >} = {aj,a^,} = 



Notations \pi.-Pn) 
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will be used. The order spin operator <r(x) 



can be determined in the infinite line x e R as a normally 



2 



ordered exponent [T3L fl4| : 



(7(X) = 

P(x) 



=P(x)/2 



dx'( X (x',y)d yX (x',y)) 



y=0' 



(4a) 
(4b) 



x(x,y) 



dp 
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where a = m 1 / 8 2 1 / 12 e _1 / 8 A 3 / 2 is the zero-field vacuum 
expectation value of the order field (spontaneous mag- 
netization), A — 1.28243... is Glaisher's constant. It is 
useful to rewrite operator : p(0)/2 : in the form 



dpdk 



- JJ-os 8TT 2 y/ui(p) u(k) 

(a p a k + a) p a\) 



(5) 



U u(k)-u(p) f uj(k)+uj(p) 



k 



One can easily verify, that Eqs. Q, J5j lead to the well- 
known formfactors ja, (pi...pjv|cr(0)|p' 1 ...p' Ar ,) of the 
order field in the infinite system. Note, that the field 
theory determined by hamiltonian @-0 is Lorentz in- 
variant. 

We are concerned with the problem of the meson mass 
spectrum in the low-temperature phase IFT in the small 
field regime ha <C to 2 . The meson masses M n are de- 
termined by the eigenvalues of the full Hamiltonian J3J): 



(H-E vac ) I $„) =M n I $„), 



(0) 



where -E V ac is the ground-state energy, and | $„) is a 
one-meson state with zero momentum. In the two-quark 
approximation one considers the projection of this equa- 
tion onto the two-fermion subspace T2 of the Fock space 
T: 



V-zHVi I #„) =M n I *„), 



(7) 



where V2 is the orthogonal projector on JF 2 . The scalar 
product of Q with the bra- vector (—p,p | yields the 
Bethe-Salpctcr equation 



(2«(p) -M„)$ n (p) = cr h 



1 



P-p' 



1 



PP 



2w(p)^) 



w(p)w(p') 

Z7T 



(8) 



first obtained by Fonseca and Zamolodchikov [6J. Here 
the principal value integral is implied. The two-particle 
wave function in the rest frame in the momentum repre- 
sentation is defined via 



* n (p) = R-V 2 (-p,p\* n ), 



(9) 



where R is the system length in the x-direction, which is 
put to infinity in the final formulas. The function $ n (p) 



is odd <E> n (— p) = — and should be normalized as 



1. 



(10a) 
(10b) 



It is worthwhile to rewrite equation (JHJ in the coordi- 
nate representation: 

(2K - M n )^ n (x) = -2o7i|x| 0„(x) + 2aW7c£ n (x). (11) 

Here <p n (x) denotes the configuration-space wave function 



</>n(x) 



dp 



e lpx $n(p), 



, 2tt 

<^„(-x) = -0„(x), 



(12a) 
(12b) 



and the integral operators K and U have the kernels 

K(x,x')= [ ^e^'V^ + m 2 , (13a) 
<ip dp' exp[i(px — p'x')] 



£/(x,x') 



4tt 2 



uj(p) - ui(p') 
p-p' 



uj(p)uj(p') 

, 1 pp' 



(13b) 



2w(p)w(p / ) 



Discrete-lattice versions of equation (|ll|l , without the last 
term in the right-hand side, were studied in papers [lil 

El. 

All terms in (|llfl have clear physical meaning. The left- 
hand side contains the relativistic kinetic energy operator 
K. The first term in the right-hand side describes the 
long-range confining attractive potential, which linearly 
increases with the distance between quarks |x|. Operator 
2ahU corresponds to the short-range diagonal interaction 
between quarks, which exponentially vanishes at large 
distances |x| 3> m. 

At small h, the lower part of the spectrum M n of equa- 
tion pit has been obtained by Fonseca and Zamolod- 
chikov 6] in the " non-relativistic" approximation 

to ~ C " 20 n + ^ V1400 280 J + 14 ' 

(14) 

where Q is a dimentionless parameter proportional to h: 
Q = 2ah/m 2 , and —z ni n — 1,2, ... are zeros of the Airy 
function, Ai(— z n ) = 0. The leading term of expansion 
11141 reproduces the old result of McCoy and Wu Q . 

It is worth stressing, that expansion (|14fl describes the 
spectrum of equation 1)11(1 , only if two conditions wee sat- 
isfied together: £ <C 1 and n < Exploiting ifHfjl and 
the large-n expansion for zeros of the Airy function 



48 



36 



+ 0(6" 9 ) 



(15) 
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with b„ — 



:>,tv 



(An - 1) 
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, one can see, that inequality 



n < C^ 1 is equivalent to QJ. Condition guarantees, 
that the wave function </> ra (x) slowly varies on the scale 
Ax ~ mT 1 . This allowed Fonseca and Zamolodchikov to 
replace the relativistic kinetic energy operator K by its 
non-relativistic counterpart, and to reduce integral equa- 
tion Ulljl to the (differential) perturbed Airy equation. 

We obtain the mass spectrum of equation (| 1 1 1) in a 
different, large-n region, namely at C *C 1, n ^> 1. The 
resulting expansion reads as 



M n =mu n + M^ +0(C 3 ), 



where u n solves equation 

1/2 



2 \ A 



arccosh(u n /2) = 



(16) 



n-^nC, (17) 



and 



( 2 m 



40 



ttd -4 



(18) 



Two notes are in order. 

(i) While the small-n expansion (|14|l is not analytical in 
£ at fixed n, the large-n expansion (f 1 Cj|> - (f 1 H|) is analytical 
in C, but at fixed n£ = const. 

(ii) The small-n and large-n regions cross at 1 <c n <c 
C _1 , covering together all possible values of n G N. In 
the crossover region 1 <C n <C C -1 , t ne both small-n and 
the large-n expansions l|14l) and l|16 fl -l|18 |) can be reduced 
to the same form: 



M„ - 2m 



-2/3 
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4/3 



20 
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-o(b-±) 



+C 



6/3 



l\bl 901 



1400 4480 



-0(C 8/3 ). (19) 
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The leading term in expansions (|16f) describes semiclas- 
sical energy levels, given by the Bohr-Sommerfeld quan- 
tization rule, as one should expect for n 3> 1. In our 
problem, the latter can be reformulated in a relativistic 
covariant form. 

Consider two interacting particles with coordinates 
Xi(t), x 2 (t) 6 K described by the classical action 

A c i = -m f" 1 di[(l-x 2 ) 1 / 2 + (l-x 2 ) 1 / 2 ]-2/iCT5', (20a) 
Jo 

with S=[ dt[x 2 (t) -xi(t)], (20b) 
Jo 

under the following constrains : xi(0) = X2(0), xi(t m ) = 
X2(i m ), and xi(t) < x 2 (t) for < i < t m . It is evident, 
that action A c \ is relativistic invariant, since S is just 
the area of the loop in the plane (x, t) , bounded by the 
particle world-paths x±(t), X2(i). Then, the quantization 
condition 



2/m>' = 2t| n- ~ ). /; - 1! 



(21) 



leads to the semiclassical mass spectrum M n = mu n of 
the pair, where u n is determined by (|17|l . Emphasize, 
that the quantization condition (|21l) is relativistic covari- 
ant, and does not imply any constrains on the total mo- 
mentum of two particles. 

Let us briefly explain now, how the large-n expansion 
Qlfi jl -i[18 jl was obtained. One could not treat the right- 
hand side of equation l|ll|) as a perturbation at a when- 
ever small h, since the product hx in the first term could 



be large at large enough x. To overcome this difficulty, 
the strategy familiar from the theory of the magnetic 
breakdown in normal metals was applied. We di- 
vided the x-axis into three intervals by the points ±a, 
where mT 1 <C a <C (rnC) 1 ; and considered (|ll|l in these 
regions separately. It is natural to call the middle interval 
—a < x < a by the "scattering region", and surrounding 
ones - by "transport regions". In the transport regions, 
the relative motion of quarks is semiclassical for n»l. 
In the right transport region x > a, the wave function 
4> n (x) approaches with exponential accuracy to the func- 
tion C(j) r (x, M n ), 



r (x,M„ 



dp 
~2~i 



exp 



i(f(p)-M n p) 



2ha 



ipx 



,(22) 



with f(p)=2 / dp' lo(p'), 
Jo 



and some constant C. One can easily verify, that func- 
tion (|22|) solves equation 

(2K- M n )(f> r {x,M n ) = -2ahx<j) r (x,M n ) (23) 

in M. In the left transport region x < —a, function 4> n (x) 
should approach to —<j) r (—x) due to the antisymmetric 
condition l|12b|l . In the scattering region —a < x < a, 
equation (|llfl was solved perturbativly in h to the first 
order. Joining at x ~ ±a solutions of different regions 
provides us the condition, which determines the discrete 
mass spectrum (ITH|l - (tTHll . 

Note, that in the zero order in h, the normalized (ac- 
cording to (|10afl . (|12afl ) wave function can be written as 



4 



<j> n (x) = J — [6(x)<fi r (x,mUn) - 9(-x)(j) r (-x, mu n )], 
V Po 

(24) 

where 9(x) is the step-function, and pq — [(M n /2) 2 — 
m 2 ] 1 / 2 . In the scattering region |x| < to -1 , the wave- 
function 1)24(1 has the asymptotics: 

0„(x) ~ (-1)" sin(pox). (25) 

Po 

The above analysis has been limited by the two-quark 
approximation J7J). However, the exact meson state con- 
tains, of course, multi-quark components, which should 
also contribute to the meson mass _6jj. Let us discuss 
qualitatively such multi-quark effects. Again, it is rea- 
sonable to consider the transport and the scattering re- 
gions separately. 

In the transport regions |x| 3> to -1 two quarks move 
independently one from another (apart from the long- 
range attractive force). Multi-quark effects here can be 
absorbed, at least partly, by radiative corrections to the 
parameters to and a in equation Ijlll) . Correction 5m 
to the quark mass known due to Fonseca and Zamolod- 
chikov 11] is of order 5m ~ ft 2 . Simple perturbative 
estimates show, that the leading correction to a is of the 
same /i 2 -order. 



In the scattering region |x| < to -1 , two quarks with 
momenta ±po collide. This collision is elastic, if M n < 
4m. In this case virtual multi-quark processes lead to the 
second-order corrections to the scattering phase, which 
in turn give rise to the meson mass correction in the 
h. 3 -order. Above the two- meson threshold, nonelastic 
scattering channels open, causing decay of mesons with 
M n > Am. The meson decay width T n can be obtained 
in the leading ft 3 -order from the Fermi's golden rule: 

T n = 2TrJ2\ (° ut I V | *„) | 2 5{mu n - E out ). (26) 

out 

Here the interaction operator V is given by ©, an out- 
state can be taken as | out) =| pi .-.p2j}, and | $„) is 
the meson state, normalized according to (|10b|> . in the 
two-quark semiclassical approximation. In equation i|26|) 
one can substitute 

since quarks collide in the scattering region, where 
asymptotics 1(251) is applicable. Then we immediately ob- 
tain from (|26|l an explicit representation for T n in terms 
of the well-known formfactors [|| 0] of the order opera- 
tor cr(0): 



r„ = 



ft 3 en 



3m 



m{u 2 



4) ^ 



dpi...dp 2j 

(27r) 2 J- 2 



5{pi-\ \-p2j)5(ui-\ \-u> 2 j-mu n ) | (P2 3 --Pi | cr(0) \po,-Pa) | 2 , (28) 



where loi = to(pi), and j m — [u n /2] stays for the integer 
part of Mn/2. The extra /i-factor appears in J2HI from 
B27JI . It reflects, that quarks spend almost all the time 
in the transport regions, and only rarely come up in the 
scattering region, from which the decay then occurs. 

In summary, we have obtained the large-n excitation 
spectrum in the Ising field theory in the weak confine- 
ment regime perturbativly in applied field h. We argue, 
that the two particle approximation is sufficient to de- 
termine the entire spectrum up to the first order in h, 
whereas, starting from the second order in h, one should 
take into account nondiagonal multi-fermion contribu- 
tions. For energies above the stability threshold, the ex- 
citation decay width has been determined in the leading 
order in h. 
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